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ABSTRACT 

We present an analysis of observations made with the Arcminute Microkelvin Imager (AMI) 
and the Canada-France-Hawaii Telescope (CFHT) of six galaxy clusters in a redshift range of 
0.16-0.41. The cluster gas is modelled using the Sunyaev-Zel'dovich (SZ) data provided by 
AMI, while the total mass is modelled using the lensing data from the CFHT. In this paper, 
we: i) find very good agreement between SZ measurements (assuming large-scale virialisation 
and a gas-fraction prior) and lensing measurements of the total cluster masses out to r2oo; ii) 
perform the first multiple-component weak-lensing analysis of Al 15; iii) confirm the unusual 
separation between the gas and mass components in A1914; iv) jointly analyse the SZ and 
lensing data for the relaxed cluster A611, confirming our use of a simulation-derived mass- 
temperature relation for parameterizing measurements of the SZ effect. 

Key words: 

cosmology: observations - cosmic microwave background - radiation mechanisms: non- 
thermal galaxiesxlusters - Sunyaev-Zel'dovich - weak gravitational lensing - galax- 
ies:clusters:individual (Abell 115, Abell 611, Abell 851, Abell 1914, Abell 2111, and 
Abell 2259) 



1 INTRODUCTION 

We employ two independent methods in a pilot study to investigate 
mass distributions of six galaxy clusters selected to cover a range 
of redshifts and merging states: 



* We request that any reference to this paper cites "AMI Consortium: 

Hurley-Walker et al. 20 11" 

t Issuing author; E-mail: nhw@icrar.org 



(i) Weak gravitational lensing (see e.g. Bartelmann & Schneider 
2001 for a review), in which images of background objects are dis- 
torted by a mass lying along the line-of-sight, and can be used to 
directly probe the cluster mass distribution. 

(ii) The Sunyaev-Zel'dovich effect (SZ; see e.g. Birkinshaw 
1999 and Caiistrom et al. 2002 for reviews), the inverse-Compton 
scattering of the Cosmic Microwave Background (CMB) by the hot 
cluster gas, which effectively measures the gas pressure. The dark 
matter content is normally assessed by combining the SZ informa- 
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tion with an X-ray measurement of gas temperature and using the 
assumption of hydrostatic equilibrium; but neither the X-ray tem- 
perature nor this assumption are necessary if we assume the cluster 
is virialised and we incorporate a sensible prior on the ratio of gas 
to dark matter (AMI Consortium: Olamaie et al. 2010). 

SZ and weak lensing have a natural complementarity, as they both 
have the potential to measure the distributed outskirts of clusters, 
with no strong bias towards concentrations of gas or mass. Also 
possible is the determination of the gas fraction of a given cluster by 
calculating the gas mass from the SZ and the total mass from lens- 
ing. It is expected that thousands of galaxy clusters will be detected 
by new SZ surveys performed by the Atacama Cosmology Tele- 
scope (ACT: Swetz et al. 2011), the South Pole Telescope (SPT: 
Vanderlinde et al. 2010) and Planck (Planck Collaboration 201 lb); 
the last has produced the first all-sky SZ catalogue, the first release 
of which is available in Planck Collaboration (2011a). Upcoming 
large-area multi-wavelength optical surveys such as the the Dark 
Energy Survey (The Dark Energy Survey Collaboration 2005) and 
eventually the Large Synoptic Survey Telescope survey (Sweeney 
2006) will improve photometric redshift measurements of galaxies 
and release lensing-quality data over tens of thousands of square 
degrees of sky. Combining SZ and lensing measurements of very 
large samples of galaxy clusters may allow us to model their inter- 
nal physics well enough for cosmological applications (Hoekstra & 
Jain 2008). 

The Arcminute Microkelvin Imager (AMI) is a radio interfer- 
ometer that has made observations of hundreds of known galaxy 
clusters to measure their gas masses and structures via the SZ. A 
limitation of the large-area SZ survey instruments is their inability 
to resolve the morphology of cluster gas, so with these instruments 
it is difficult to examine the cluster dynamical state. Interferomet- 
ric arrays such as AMI and the Combined Array for Research in 
Millimeter-wave Astronomy (Bock et al. 2006) allow the examina- 
tion of structures on scales between the high resolution of X-ray 
instruments and the somewhat lower resolution of the SZ survey 
instruments, as well as covering the northern sky. 

This study concerns a small selection of clusters of known 
X-ray structure, observed by AMI, for which publicly-available 
CFHT data were accessible at the time the AMI observations were 
made. As a pilot study for a future SZ-lensing comparison with a 
larger sample, we examine here six clusters with: a redshift range 
of 0.16-0.41, arange of masses and varying degrees of merging ac- 
tivity. Given the depth of the optical observations (see Section 3), 
clusters at higher redshift would be more difficult to observe as 
the field galaxy selection would likely be contaminated with fore- 
ground galaxies, and with only two optical bands we would be un- 
able to easily reduce this contamination. Below a redshift of about 
0.08, AMI starts to resolve out the cluster gas; this sample of clus- 
ters should not be affected by either of these issues, allowing us 
to examine the agreement of lensing and SZ mass measurements 
and the effect of the cluster dynamical states. AMI is limited to 
observing above Declination 20° and is effectively limited to mea- 
suring clusters of internal gas temperatures ^ 2keV. Clusters of 
lower temperature may be detectable but would require very high 
(unphysical) electron densities to produce detectable SZ flux. 

A data reduction pipeline was developed to extract a weak- 
lensing catalogue from background galaxies in the CFHT fields. 
This allowed the measurement of the total matter distributions of 
the galaxy clusters, which could then be directly compared with 
gas measurements from the SZ. The SZ observations and data re- 
duction are described in Section 2, and weak lensing in Section 3. 



The Bayesian analysis, including sampling parameters and priors, 
is described in Section 4. We present notes on each cluster and its 
available data in Section 5, discuss the ramifications of these results 
in Section 6 and outline our conclusions in Section 7. 

Throughout, we assume a 'concordance' ACDM cosmology, 
with £2 m o = 0.3, £2 A ,o = 0.7, and Hq = 70kms~ 1 Mpc~ 1 (and thus 
h = Hq/ 100 = 0.7). All coordinates are J2000 epoch and all optical 
magnitudes use the AB system. 



2 SZ OBSERVATIONS 

AMI is a dual set of aperture-synthesis arrays located at the 
Mullard Radio Astronomy Observatory, Lord's Bridge, Cambridge, 
UK. The AMI Small Array (SA) consists of ten 3.7-m-diameter 
equatorially-mounted dishes with a baseline range of ~ 5-20 m, 
while the AMI Large Array (LA) has eight 12.8-m-diameter dishes 
with a baseline range of ~ 20-100 m. Both arrays observe I + Q po- 
larisation flux densities in the band 12-18 GHz, each with system 
temperatures of about 25 K. 

The back ends are analogue Fourier transform spectrometers, 
from which the complex signals in each of eight channels of 750- 
MHz bandwidth are synthesised, and the signals in the synthesised 
channels are correlated at the ~ 10 per cent level. In practice, the 
two lowest-frequency channels are generally not used due to in- 
terference and a poor correlator response in this frequency range. 
Further details are given in AMI Consortium: Zwart et al. (2008). 

Details of the observations for these six clusters, made in 
2008-2009, are shown in Table 1. The SA observed a single 
pointing while the LA used a (61+19)-point raster mode with 
4 arcmin spacing. Phase calibrators were chosen from the Jodrell 
Bank VLA Survey (Patnaik et al. 1992; Browne et al. 1998; Wilkin- 
son et al. 1998) on the basis of proximity (^ 5°) and flux density 
(> 1 Jy at 15 GHz). 

The AMI data reduction was performed using our in-house 
reduction software REDUCE. This is used to apply path-delay cor- 
rections, to flag interference, shadowing and hardware errors, to ap- 
ply phase and amplitude calibrations and to Fourier transform the 
correlated data to synthesise the frequency channels, before writ- 
ing output wv-FITS files suitable for imaging in the Astronomical 
Image Processing System (MPS; Astronomical Image Processing 
Software 2007). 

Flux calibration was performed using short observations of 
3C48 and 3C286 near the beginning and end of each run, with as- 
sumed I + Q flux densities for these sources in the AMI channels 
consistent with Baars et al. (1977) (see Table 2). As Baars et al. 
measure I and AMI measures I + Q, these flux densities include cor- 
rections for the polarization of the sources. After phase calibration, 
the phase of AMI over one hour is generally stable to 5° for chan- 
nels 4—7, and to 10° for channels 3 and 8. The system temperatures 
of each AMI antenna are continously monitored using a modulated 
noise signal injected at each antenna; this is used to continuously 
correct the amplitude scale in a frequency-independent way. The 
overall consistency of the flux-density scale is estimated to be bet- 
ter than five per cent. 

To transform the SA data into a format suitable for our anal- 
ysis, the unsmoothed wv-data for all good observations were con- 
catenated together to make a single wv-FITS file for each channel. 
These were then transformed into lists of visibilities for the purpose 
of generating a covariance matrix for the data. 

The FITS data were also binned into cells of width 40A. This 
was determined by Marshall et al. (2003) as a suitable binning scale 
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Figure 1. LA raster signal-to-noise maps of A115, A1914, A2111, A2259, A611 and A851, primary-beam-corrected, with sources detected at 4o~ and above 
labelled with identifying letters which refer to the sources as shown in Tables 7-12. Greyscale and contours are both signal-to-noise, with the latter in levels 
of +3, +4, +5, +6, +8, +10, +20, +30, +50 and +100. 
© 2011 RAS, MNRAS 000, 1-21 
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Cluster RA Dec z S A observing SA map noise/ LA observing LA 19-pt map noise / LA 61 -pt map noise / 



(J2000) (J2000) time/hrs ^iJybeam 1 time/hrs jiJybeam 1 jiJybeam 

ATl5 00 55 59.50 +26 19 14.0 (1197 67 114 19 101 23T~ 

A1914 14 26 02.15 +37 50 05.8 0.171 47 115 22 111 291 

A2111 15 39 44.08 +34 24 56.2 0.229 43 90 22 99 268 

A2259 17 20 10.60 +27 40 08.4 0.164 52 100 18 113 248 

A611 08 00 59.40 +36 03 01.0 0.288 76 85 23 76 208 

A851 09 43 07.08 +46 59 51.0 0.410 77 75 28 70 188 



Table 1. Details of the AMI observations of the six clusters. 



Channel v/GHz S(3C286)/Jy 5(3C48)/Jy FWHM LA / arcmin 
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13.9 


3.74 


1.89 6.08 


4 


14.6 


3.60 


1.78 5.89 


5 


15.3 


3.47 


1.68 5.70 


6 


16.1 


3.35 


1.60 5.53 


7 


16.9 


3.24 


1.52 5.39 


8 


17.6 


3.14 


1.45 5.25 



Telescope (CFHT). MegaCam offers a field-of-view of one square 
degree, which is more than sufficient for examining the lensing sig- 
nal from our six clusters. After retrieval of the archival data they 
were astrometrically and photometrically calibrated and finally 
co-added as described in Erben et al. (2009), using textscswarp 
(seehttp: //www. astromatic .net/ software/ swarp). This re- 
sulted in approximately degree-square frames with pixel sizes of 
0.186 arcsec. 

We give important characteristics of the finally co-added data 
in Table 3. The limiting magnitude in this table is defined as 
the 5-ct detection limit in a 2 arcsec aperture via m iim = ZP — 
2.51og(5^//Vpi^a s ]jy), where ZP is the magnitude zeropoint, N p [ x 
is the number of pixels in a circle with radius 2 arcsec and rj^y the 
sky background noise variation. S EXTRACTOR (Bertin & Arnouts 
1996) was run using default parameters on each of the images to 
extract all sources. Generally for each image, around 10 s objects 
were detected in the r-band, and around 7 x 10 4 in the g-band. An 
object was defined as consisting of 10 or more pixels above 3 x the 
local noise. 



3.1 Point spread function (PSF) 

The point spread function (PSF) describes the convolution of the 
image due to the blurring effect of the atmosphere, telescope wind 
shake, telescope optical distortions and CCD diffusion. It is essen- 
tial to correct for the effects of the PSF, which can be at least as 
strong as the lensing signal. In these CFHT images, the PSF ellip- 
ticity distorition is of the order of a few per cent, and is not always 
restricted to the edges of the field. Usefully, foreground stars are 
unresolved and trace out the PSF, so by measuring stellar elliptici- 
ties and sizes we can measure and then correct for the PSF. 

Optically-detected objects can be described by five param- 
eters, their positions x, y, the ellipticity vectors £i, £2 and ab, 
the product of the semi-major (a) and semi-minor (b) axes at the 
FWHM and thus describing the size of the objects. The average 
FWHM of the stars for an image is taken as the 'seeing', and also 
has an effect on the lensing signal: poor seeing reduces the ellip- 
ticity of objects and makes them appear more round: this tends to 
circularise background galaxies and thus weaken the shear signal. 
£1 and £2 are the spin-2 tensors: £1 = ecos(20); £2 = esin(20), 
where e is the ellipticity (a — b)/(a + b) and 6 is the orientation, 
measured clockwise from south. 

Within the source catalogues, stars form distinct populations 
of low-FWHM objects and were easily extracted. Stars with SEx- 
TRACTOR flags of less than 4 (i.e. their shapes are well-known and 
there are no nearby contaminating objects) were used as the tem- 
plate stars for fitting the PSF. Typically for these images, around 
3 x 10 3 stars were visible, with the exception of A2259, whose lo- 
cation near the Galactic plane resulted in the detection of 10 4 stars. 
The r-band images of these stars were further analysed using the 



Table 2. Assumed I + Q flux densities of 3C286 and 3C48 over the 
commonly-used AMI band, and the AMI LA primary beam Full-Width- 
Half-Maxima (FWHM) for each channel. 

which reduces the size of the data to a manageable level without 
adversely affecting the resulting inference of cluster properties. 

Maps were made using imagr in AIPS for each channel of 
the SA and LA; however we present only the combined-channel 
maps of the SA and LA observations, imean was used on the LA 
individual maps to attach the map noise to the map header, and 
f latn was then used to stitch the maps together, with a primary 
beam correction applied using parameters shown in Table 2. 

Source-finding was carried out on the LA data using software 
modified from that used in the 9C survey (Waldram et al. 2003). 
Spectral indices were fitted using LA maps for all six channels, as- 
suming source fluxes, S, follow a power-law relation of 5 V~ a for 
the AMI frequencies, v. The properties of point sources detected at 
4<r and above by the LA are used as priors when modelling the 
SA data in the analysis (see Section 4). Fainter sources more than 
5 arcmin from the pointing centre are directly subtracted in order to 
reduce the computational time needed; these are marked with a ' x ' 
instead of a '+' in the AMI maps. 

The position prior is a delta function since the resolution of 
the LA is around three times that of the SA. The flux densities S 
are given Gaussian priors: a is given by a conservative calibration 
error of five per cent added in quadrature to the local map noise. 
Spectral index (a) priors are also Gaussians with a equal to the 
error on the spectral index fit. These errors tend to be small (~ 0.4) 
for bright sources and large (~ 2) for faint sources. 

Fig. 1 shows LA maps for the six clusters with the detected 
radio point sources labelled with identifying letters. The positions 
and properties of the sources are given in Tables 7-12, sorted by 
measured LA flux density. 



3 WEAK GRAVITATIONAL LENSING DATA 

The clusters observed by AMI are of a few arcminutes angular ex- 
tent. High-quality optical data for the six clusters were retrieved 
from the public Canadian Astronomy Data Centre archive (see 
http://cadcwww.dao.nrc.ca/cadc/). All clusters were ob- 
served for weak-lensing in the r-band (A e ff = 623 nm) and for sup- 
plementary colour information in the g-band (A e ff = 477 nm), us- 
ing MegaCam (Boulade et al. 2003) on the Canada-France-Hawaii- 
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Cluster 


obs. dates 




P.I. 


filter 


exp. time [s] 


seeing [arcsec] 


ZP mag. [AB 


A115 


08/2004 - 10/2004 


H. 


Hoekstra 


r 


6602.07 


0.69 


24.87 


A115 


08/2004 - 10/2004 


H. 


Hoekstra 


8 


1600.65 


0.79 


23.53 


A1914 


05/2006 


H. 


Hoekstra 


r 


6001.9 


0.71 


24.56 


A1914 


05/2006 


H. 


Hoekstra 


8 


3601.63 


0.90 


24.97 


A2111 


05/2006 - 06/2006 


H. 


Hoekstra 


r 


6482.13 


0.66 


24.94 


A2111 


05/2006 - 06/2006 


H. 


Hoekstra 


8 


1800.75 


0.55 


24.34 


A2259 


08/2004 


H. 


Hoekstra 


r 


4001.29 


0.82 


24.18 


A2259 


08/2004 


H. 


Hoekstra 


8 


1600.65 


0.82 


24.00 


A611 


12/2004-01/2005 


H. 


Hoekstra 


r 


4801.59 


0.77 


24.29 


A611 


12/2004-01/2005 


H. 


Hoekstra 


g 


2520.93 


0.82 


24.52 


A851 


11/2004 


H. 


Hoekstra 


r 


6602.67 


0.95 


24.40 


A851 


11/2004 


H. 


Hoekstra 


g 


2751.05 


0.77 


23.89 



Table 3. Quality information on the co-added optical MegaCam data. 



program IM2SHAPE, which fits a Gaussian to each object using a 
Bayesian method to find the most likely parameters (Bridle 2001). 

The r-band images of Al 15, A2259 and A61 1 had discontinu- 
ous behaviour along the joins between the individual frames. These 
joins, and areas directly around bright stars and their diffraction 
patterns, were excluded from the final lensing analysis, reducing 
the available image areas by ~ 20 per cent. These regions are never 
directly over cluster centres, are less than 0.5 arcmin wide and are 
in a gridlike pattern so the decrease in available field galaxies does 
not strongly affect any particular radial bin. Objects with extreme 
(e > 0.95) ellipticities were also discarded. 

The functions of £\ , £2 and ab for the stars contained noise 
components dependent on the individual frames and the SWARP 
process, so the functional forms were not known. Given the gaps 
discussed above and the inclusion of this noise component, a sim- 
ple interpolation scheme was not optimal. In order to fit a smooth 
and continuous form for each function an artificial neural network 
was used to 'learn' the noise component and ultimately remove it, 
and interpolate over regions of sparse or no data. Details of the net- 
work software used can be found from a previous application to 
cosmology in Auld et al. (2008). 

From these functions, a reference file describing the PSF was 
generated for the background galaxy positions, to be deconvolved 
when fitting the galaxy shapes. Residuals of £\ and £2 for the stellar 
catalogues after PSF fitting (calculated simply as ej tar — eJ SF ) are 
shown in Fig. 2. 



3.2 Contaminant and cluster galaxy removal 

Once PSFs had been determined for each r-band image, IM2SHAPE 
was run again, this time on the background galaxy catalogue, ex- 
tracted from the initial full set of objects by: 

• removing the stars and saturated objects, obvious from their 
FWHMs; 

• removing very blue objects, generally foreground galaxies 
with large angular extent; 

• removing the cluster galaxies; 

• making a brightness cut in the r-band: r > 21. 

Cluster red sequence galaxies were found from their overdensity 
in the (g — r):r colour-magnitude diagram and examining those ob- 
jects in the x-y plane. We observed dense spatial clustering at the 
cluster positions indicating that the cut included the cluster galax- 
ies. Their distributions on the sky also give rough estimates of the 
angular extents and thus sizes and shapes of the galaxy cluster 
mass distributions, although we do not attempt to use this informa- 



tion quantitatively in this paper. Table 4 shows the colour proper- 
ties of the galaxies lying in the red sequence region and including 
the cluster members, and Fig. 3 shows the selections in (g — r):r 
space. A check on the selection method was made by examining 
the spectroscopically-confirmed cluster members for A115 (Bar- 
rena et al. 2007) and A2111 (Miller et al. 2006 and Ascaso et al. 
2008) and matching their positions with those in our catalogue; they 
were found to lie along the same red sequences. Slightly older data 
from Dressier et al. (1999) and Belloni & Roeser (1996) were used 
to find likely cluster members of A851, with an estimated reliabil- 
ity of 80 per cent. The brightest cluster galaxies of the three less 
extensively-observed clusters A1914, A611 and A2259 are also 
marked in Fig. 3, originally identified by Stott et al. (2008). 

For each cluster, around 2 x 10 3 objects were excluded for 
being on the red sequence, and a few hundred contaminating 
large, blue, or very bright objects were removed. After a cut at 
r-magnitude= 21 to exclude the brightest objects, catalogues of 
5-8 x 10 4 'background' galaxies were therefore obtained for each 
cluster, corresponding to number densities of 14-22 galaxies per 
square arcminute. 

3.3 Background galaxies 

IM2SHAPE was run on the background galaxies, and its output used 
to generate a shear catalogue of x, y, £\ , £2, o" ei and a £l . The er- 
rors o" ei and a £l are the fitting errors of the respective ellipticity 
vectors added in quadrature with the root-mean-squared (r.m.s.) of 
the intrinsic ellipticity distribution of galaxies, 0.25 (Hoekstra et al. 
2000). 

A simple test of whether the PSF deconvolution has been per- 
formed correctly, and the lensing signal is the dominant component 
of the map, is to bin and compare the E- and B-modes of the el- 
lipticities: £, and £ x , respectively. The former should be positive 
and decreasing with angular distance from the cluster centre, and 
the latter should be zero within the errors (see e.g. Schneider et al. 
2006). This test was performed by binning the ellipticities in radial 



Cluster 


intercept at g - r 


gradient 


max /--band magnitude 


A115 


2.50 


-0.0677 


23.0 


A1914 


2.24 


-0.0573 


23.2 


A2111 


2.95 


-0.0918 


23.2 


A2259 


2.48 


-0.0768 


23.5 


A611 


3.44 


-0.1007 


23.1 


A851 


3.06 


-0.0718 


23.2 



Table 4. Colour properties of the red sequences for each galaxy cluster, 
plotted in Fig. 3. 
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(d) A2259 (Stott et al. 2008). (e) A61 1 (Stott et al. 2008). (f) A85 1 (Belloni & Roeser 1996). 

Figure 3. Red sequences identified by overdensity in the g — r:r colour-magnitude plot. Reliably spectroscopically-identified cluster galaxies are circled and 
references are given in the subfigure captions. In the case of A851, the identification of cluster members is only accurate to 80 per cent (Belloni & Roeser 
1996) and these are encased by Ds. In each subfigure, the solid line shows the red sequence fit described in Table 4, while the dashed lines show the upper and 
lower limits on the space excluded from the field galaxy selection. 



bins of width 3 arcmin, from the fitted cluster centroid positions 
(Section 5) for A21 11, A61 1 and A1914. For Al 15, the centre used 
was the position of the highest-mass clump; for A851, two separate 
mass overdensities were detected, so the signal is shown for both 
of the components. The plot for A2259 is not shown as the lens- 
ing signal was too weak to recover good cluster parameters. Fig. 4 
shows the binned ellipticity vectors and, as expected, the B-mode 
signal is negligible given the errors, calculated as the r.m.s. of each 
bin divided by the square root of the number of galaxies in that bin. 



To model the mass of the cluster via lensing, we needed 
some estimate of the redshifts of these lensed field galaxies. For- 
tunately, the CFHT Legacy Survey Deep Field catalogues (see 
http://www.cfht.hawaii.edu/Science/CFHTLS/) were pro- 
duced using the same instrument and filters as our data, but to a 
greater depth, and with much greater redshift accuracy. We were 
therefore able to make cuts within this catalogue at the same levels 
of those of our own data, and then average the photometric redshifts 
of the galaxies within this selection. This resulted in a redshift dis- 
tribution with a mean of 1.2, a median of 1.1 and an r.m.s. of 0.7. 
With the Bayesian analysis method we employ (see Section 4) we 
are able to incorporate this large uncertainty as a sampling param- 
eter, propagating the uncertainty through to the estimation of our 
cluster parameters. 



4 MODELLING AND PRIORS 

To analyse the AMI cluster observations we use a Bayesian anal- 
ysis methodology (Marshall et al. 2003 and Feroz et al. 2009b) 
which implements MultiNest (Feroz & Hobson 2008 and Feroz 
et al. 2009a), an application of nested sampling (Skilling 2004), 
to efficiently explore multi-dimensional parameter space and to 
calculate Bayesian evidence. This analysis has been applied to 
pointed observations of known clusters (e.g. AMI Consortium: 
Rodriguez-Gonzalvez et al. 2010; AMI Consortium: Zwart et al. 
2010), and also to detect previously unknown clusters (AMI Con- 
sortium: Shimwell et al. 2010). As is standard in Bayesian meth- 
ods, priors are given for sampling parameters; the LA data is used 
to produce priors for the radio sources present in the SA data. 



4.1 SZ 

The SZ effect from a cluster can be measured by its comptonisation 
parameter, y, which is the integral of the gas pressure along the line 
of sight / through the cluster: 

y(s) = / n e k B Tdl, (1) 

where Oj is the Thomson scattering cross-section, T is the temper- 
ature of the ionised cluster gas, w e is the electron number density, 
(measured using Equation 7), m e is the electron mass, c is the speed 
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distance to the cluster. A cluster with comptonisation parameter y 
appears as a surface brightness fluctuation of magnitude 



8l v =f(v)yB v (T CMB ), 



(2) 



where £ v (7cmb) ' s tne CMB blackbody spectrum and /(v) is the 
SZ spectrum. Interferometers such as AMI effectively measure spa- 
tial fluctuations in surface brightness 8l v in the visibility plane, so 
without any transformation to the map plane (and without the as- 
sociated problems wherein), the model prediction can be directly 
compared to the data. 

The likelihood function for a set of cluster parameters x can 
be written: 



L sz = Pr(d | x) = ^-exp 



-r_ 

2 



(3) 



where % 2 ls a statistic quantifying the misfit between the observed 
data d and predicted data d p (the latter of which is a function of the 
model SZ surface brightness Sl v ): 



X 2 = (d-dP) / C-'(d-dP), 



and the normalisation factor is: 

Z L = (2w) w * 



(4) 



(5) 



which describes the terms that contribute to the data but are not part 
of the model: the CMB, thermal noise from the telescope and con- 
fusion noise from unresolved point sources. The last were modelled 
using the Tenth Cambridge Radio Survey (10C: AMI Consortium: 



Franzen et al. 201 1; AMI Consortium: Davies et al. 201 1), integrat- 
ing the confusion power from zero up to the source detection limit 
on the LA rasters. 

The cluster geometry, as well as two linearly independent 
functions of its temperature and density profiles, must be spec- 
ified in order to compute the Comptonization parameter. For the 
cluster geometry we test both spherical and elliptical profiles; the 
elliptical model is not properly triaxial, merely stretching out the 
model via an ellipticity parameter 77 at an angle tj> (measured anti- 
clockwise from west). The temperature profile is assumed to be 
constant throughout the cluster and a j3 -model is assumed for the 
cluster gas density, p g (Cavaliere et al. 1978): 



PgW : 



Pg(0) 



1 + 



(if 



(6) 



where 



p g (r)=H e n e {r), 



(7) 



jl e — 1 . 14m p is the gas mass per electron and m p is the proton mass. 
The core radius r c gives the density profile a flat top at low y and 
p g has a logarithmic slope of 3/3 at large —. 

X-ray temperature measurements can be used to break the de- 
generacy in the SZ y parameter between gas mass and temperature, 
but these are often overestimates of the temperature of the gas de- 
tected via SZ, if, as is usually the case, they are made at much 
smaller radius. Instead we use a theoretical mass-temperature rela- 
tion to constrain the degeneracy: this was introduced in AMI Con- 
sortium: Rodriguez-Gonzalvez et al. (2010) and assumes that the 
cluster is virialised, and that all kinetic energy is in the internal 
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energy of the cluster gas; it does not assume hydrostatic equilib- 
rium. AMI Consortium: Olamaie et al. (2010) show that this is a 
useful model for SZ observations, testing it on simulated data and 
extracting physical cluster parameters, and finding these to be in 
agreement with true input values in simulations. Thus we have: 



, rp Gp e M 
2^200 



G/x 



(4^(200p C ri.)y 



= 8.2keV( - 

V 



M 



2/3 



W^h^Mr. 



I©, 



'H(z)\ 2/3 
.Ho J 



(8) 
(9) 

(10) 



The center of the cluster profile is also fitted for via two position pa- 
rameters, x and y; the former equal to negative Right Ascension, in 
arcseconds, and the latter equal to Declination in arcseconds. Thus 
the final gas model has seven parameters (nine in the case of ellip- 
tical geometry): x, y, (r/, <p), p\ r c , / ga s(>"20o) and M T (r 2 oo)- The 
last parameter is the total mass inside a radius, r2oo, at which the 
total density is 200 x p cr j t (z), the critical density for closure of the 
Universe at redshift z. For simplicity, this mass will henceforth be 
denoted M, and it was given a log uniform prior varying between 
10 13 - 2 x 1O 15 /! _1 M , which is a physically reasonable range of 
cluster masses. Gaussian priors were used for the position param- 
eters, centred on the cluster catalogue positions in Table 1 with a 
set to 1 arcmin. J], (j>, j8 and r c were given uniform priors between 
0.5-1, 0-180°, 0.3-2.5, and 10-1000 /r'kpc, respectively. 

The prior on /gasC^Oo) was set to a narrow Gaussian centred 
at the 90 per cent of the Wilkinson Microwave Anisotropy Probe 
7-year (Komatsu et al. 2011) best-fit value of the baryonic mass 
fraction, / g (r 2 oo) = 0.123/r 1 , with rj = 0.02k- 1 . 

Sample visibilities are generated for the model within the prior 
ranges, compared to the data and the process is iterated until the 
sampling converges and no further improvement can be made. 



4.2 Weak gravitational lensing 

The model used for dark matter distributions throughout this analy- 
sis is the Navarro-Frenk- White profile (NFW: Navarro et al. 1995): 



PmassM _ &c 
Pcrit(z) ~ (r/r s ) + (l + r/r s ) 2 



(11) 



where r s is a scale radius and S c is the characteristic overdensity 
of the halo and is related to the concentration parameter c. The 
equation can be usefully rewritten in terms of a parameter x related 
to the virial radius r2oo by x = r/r2<x>: 



Peril (z) CX+ (1 +cx) 2 ' 



(12) 



The lens potential y for the NFW profile is derived in Meneghetti 
et al. (2003): 



w- 



4p s f" s 



, u n-d 

2 arctanhi/ 

I V l+d 



(13) 



where d is the projected radius scaled by the scale radius: d = 
9/(r s /Dg), £ cr jt is the critical surface mass density for strong lens- 
ing and p s is a scale density. For an elliptical NFW profile, d is 



scaled by the ellipticity e as 1 : 



d^td' = \ d 2 {l-e) + - 



dl 



l-e) 



(14) 



In a similar way to the use of the Comptonisation parameter y in the 
SZ analysis, the convergence K and shear y are generated at each 
galaxy position 0! , 82 using the lens potential: 



1 d 2 W d 2 V 

2 { de 2 + de 2 



7i 



12- 



d 2 \\l d 2 \ff 



de 2 



do 2 



d 2 w 



(15) 



(16) 



(17) 



and compared to that generated by a model set of cluster parameters 
x in a similar way to the SZ likelihood (Equation 3) 



Pr(d|x) 



1 

Z~L 



exp 



r 
2 



(18) 



where y} is the misfit statistic testing the measured lensed elliptic- 
ity components £,-, taken as having been drawn independently from 
a Gaussian distribution of N lensed galaxies with mean gj and vari- 
ance cr: 



i=\j=\ 



The normalisation factor Zi is 



2ko< 



The mass contained within the scaled radius d is 



M(d)=4 W 3 



log(l+d) 



l+d 



(19) 



(20) 



(21) 



and when r = 7-200, d becomes c, the concentration parameter. More 
detailed derivations are given in Marshall et al. (2003). 

Thus the mass model we use has five parameters (seven in the 
case of an elliptical profile): x, y, (77, 0), M, c, and Zfi e id> me l ast °f 
which is incorporated into Dg . These model statistics are immedi- 
ately useful when comparing to data from the literature and our SZ 
measurements. The geometry and M priors were the same as for the 
SZ modelling; c was given a uniform prior between 0.1-15, cover- 
ing a physically reasonable range. Zfr e \d was given a Gaussian prior 
centred on 1.2 with o" =0.7, as discussed in Section 3. Our data 
put little constraint on this parameter so the posterior resembles the 
prior, and thus is not plotted in the results (Section 5). 

Fitting multiple components is performed for disturbed clus- 
ters; the highest-evidence model is discussed in the text. For A61 1, 
where the large-scale gas and dark-matter distributions were re- 
laxed, it was also possible to run a joint analysis in which the mass 
was constrained by the lensing data, and the gas fraction was given 
a uniform prior. 



1 We note that the published version of Meneghetti et al. (2003) contains a 
version of this equation that is in error, and show the corrected version here. 
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5 RESULTS 

Full unmarginalised posteriors are available from the analysis; 
however, these are multi-dimensional and too large to display. In 
order to highlight particular features, small sets of parameters are 
shown plotted against each other: for instance, axis ratio Tj against 
angle <j>. The contours contain posterior probabilities at 68% and 
99%. The mean value for each parameter is given in tables in the 
section discussing each cluster, with errors at the 68% level. 

Source-subtracted SA maps are generated by using the AIPS 
task uvsub on the wv-data. The source parameters used are the 
mean source fluxes and spectral indices generated by fitting source 
models to S A data. LA data are only used to provide priors for these 
fits. The source-subtracted maps are solely for display purposes; all 
cluster parameter fitting is performed in the wv-plane. 

A useful tool for displaying the likely mass distributions caus- 
ing shear in lensing data is LensEnt (Marshall et al. 2002). While 
it can also be used as a tool to extract cluster parameters, it is used 
here only to reconstruct the mass distributions for display purposes. 
We now examine the SZ and lensing measurements for each cluster 
in detail. 



5.1 A115 

Al 15 is the most disturbed cluster in the sample: it is fully bimodal 
in X-rays, with a brighter, triangular clump in the north and a dim- 
mer, elliptical clump in the south. The two distinct gas clumps are 
distorted and in motion both tangentially and along the line-of- 
sight (Gutierrez & Krawczynski 2005). Govoni et al. (2001) also 
note that this cluster has a bright relic radio source trailing from 
the northern clump out to the north-east of the cluster. The source 
coincident with the highest X-ray emission of the northern clump 
is identified as source D in our data (Table 7; Fig. 1). At higher 
resolution, this source has a double radio structure and has been 
studied in detail by e.g. Giovannini et al. (1987). From NVSS data 
at 1.4 GHz, we find an integrated flux density of 1.35 ± 0.07 Jy, 
whereas at 15.7 GHz it has fallen to just 5.4 ± 0.3 mJy, implying a 
spectral index of a = 2.28 ± 0.03, in agreement with a = 2.5 ± 0.4 
from the LA separate-channel data alone. 

The relic radio structure is resolved out by the LA, but one can 
see immediately from the LA point-source map (Fig. 1) that the re- 
gion contains many radio point sources, including two very bright 
sources of 55 and 26mJy (sources A and B). Attempts were made 
to model the cluster and sources;unfortunately, the region also con- 
tains ~ 25 sources, leading to a very high-dimensionality problem. 
The high density of sources on top of the X-ray gas also means that 
there is a large probability space to explore when modelling the 
source and SZ flux, so the algorithm is very slow to converge on 
a solution. One problem that was immediately noticeable in early 
attempts was that source B appeared to have a much higher flux 
density in the SA data than in the LA data. This is unusual because 
in the other cluster analyses, the agreement between LA fluxes and 
SA fluxes is very good. 

Looking at the SA map of Al 15 (Fig. 5a), source B is boosted 
to 36mJy after primary-beam correction, compared to 27.6 mJy 
measured by the LA. This could be for two reasons: a) dimmer 
sources in the surrounding region are unresolved so their flux den- 
sity contributes to the flux density of source B in the map plane; 
b) source B has extended structure that is resolved out by the LA, 
resulting in a lower apparent flux density. If the first point were a 
problem, the flux density allocated to the other surrounding sources 
would be increased in order to fit the data. However, in the uv- 



plane, the power does not appear to be coming from these sources, 
as there is no evidence for these sources having flux densities larger 
than those measured by the LA. Only source B appears to have a 
boosted flux density. This implies that the source is extended, which 
is not surprising given the extended relic structure visible at lower 
frequencies. There was no variation in the flux densities measured 
from observations at different dates, so the difference is not due 
to variability. The conclusion is that the source is an ellipsoid ex- 
tended roughly north-south, with an extent larger than the SA syn- 
thesised beam. Of course, extended positive sources are degenerate 
with the SZ signal of the cluster gas, so it would be unlikely that 
an accurate SZ signal could be extracted from the sky immediately 
behind and around this source. 

A negative signal is visible around the location of the south- 
ern X-ray clump (00 h 55 m 59?5, +26°20'02") in the SA source- 
unsubtracted data with a peak negative flux of — 1.4 mJy (~ lOo"). 
Some of the decrement may be created by improperly CLEANed 
sidelobes of nearby sources, particularly source B. However one 
does not see negative signal on all sides of the positive source con- 
glomeration, so it is likely that a large component of this negative 
signal is SZ. Unfortunately, the data are insufficient to model the 
gas content of the cluster. 

The LensEnt mass reconstruction is shown in Fig. 5b and is 
in visual agreement with that made by Okabe et al. (2010a), with 
five mass clumps resolved; their positions are overlaid. We also fit- 
ted elliptical NFW profiles to the data; the model with the highest 
evidence contained five components and the resulting parameters 
are given in Table 5. The position posteriors of the fits to the five 
profiles were examined for evidence that they were composed of 
yet further components, but it appears that the position posterior 
distributions, visible in Fig. 5c, are unimodal. Okabe et al. mea- 
sure M200 by fitting a NFW profile to the central mass peak and 
found M 2 oo = 4.45+j 35 x 1O 14 /i _1 M , the upper limit of which is 
in agreement with our measurement of component 3, in which the 
bulk of the mass is found, and which is coincident with the southern 
X-ray clump. This mass is highly unconcentrated, with c < 1; this 
is in visual agreement with the LensEnt reconstruction, but may 
show that the NFW profile is not well-suited to fitting this irregular 
distribution. 

Interestingly, component 1 is has a similar mass to compo- 
nent 4, despite only the latter clump coinciding with X-ray emis- 
sion. This implies that the bulk of the X-ray emission in the north- 
ern clump is from the X-ray-emitting double radio source, and there 
may not be a large dark matter component to this clump. Com- 
ponent 1 might be a bullet-cluster-like dark-matter blob that has 
passed through the cluster leaving behind its stripped gas. An el- 
liptical galaxy at 00 h 56 m 19?7, 26°21'53" is within the mass over- 
density of component 1. It lies on the red sequence of A 1 15, which 
spatially overlaps the area of component 1. 

Component 2 is clearly detected and it is uncertain what rela- 
tion this object bears to the bulk of the cluster. At a total mass of just 
~ 9 x 10 13 /z _1 M Q , its gas mass should be of order ~ 1O 13 /! _1 M , 
rendering it undetectable by AMI in the SZ. Visual examination of 
the r-band image shows no unusual galaxy overdensity at this po- 
sition, and galaxies there do not seem to lie on the red sequence of 
A115. Spectroscopic follow-up, or observations in another optical 
band to compare the densest region of galaxies of A115 with this 
area might indicate whether this mass concentration is sub-cluster- 
like, related to A115, or an unrelated mass overdensity along the 
line-of-sight. Component 5's position is not well-constrained com- 
pared to the other components and it has a low concentration and 
uncertain mass. 
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Figure 5. A115. The SA map is shown without source-subtraction; contours are every 3a, and truncate at +30c for clarity. The parameters of the labelled 
sources can be found in Table 7. In this and all subsequent SA and LA maps in this section, the CLEAN process has been applied down to 3o~, no correction 
for primary beam has been applied, and the synthesised beam is represented by an ellipse in a box in the lower left-hand corner. The second image shows 
a LENSENT 2.5arcmin reconstruction of A115 overlaid with the positions and labels of the components (Component x =Cx) whose fitted parameters are 
given in Table 5. The greyscale is linear and is overlaid with 20% density contours. The relative lengths of the cross arms indicate the axis ratios 7], and the 
inclinations indicate the angles <j> of the components. The final image shows the A115 weak-lensing 2D posteriors; common ranges of x,y = ±500 arcsec, 
c = 1 - 15, and M = 10 13 - IO^/t'Mq are shown for all components except component 3, and delineate the prior space searched for the five components 
found. The posteriors for rj and <p are not shown as they are only well-constrained for component 3. 



Component 


1 


2 


3 


4 


5 


x 1 arcsec 


-327 ±6 


-200 ±10 


-49 ±9 


161 ±8 


355 ±44 


y 1 arcsec 


169 ±5 


-328 ±8 


133 ±24 


389 ±9 


-394 ±51 


<t>/° 


85 ±60 


100 ±50 


123 ±3 


75 ±53 


1 1 3 ± 45 


V 


0.88 ±0.08 


0.85 ±0.09 


0.71 ±0.07 


0.87 ±0.07 


0.83 ±0.09 


c 


8.8± 1.7 


5.3±0.9 


0.4 ±0.1 


8.0±1.0 


9.8 ±2.5 


M IU i; /; 'M. 


0.7±0.2 


0.9±0.2 


9.3 ±0.7 


1.3±0.2 


0.4 ±0.1 



Table 5. Mean posterior values for the five elliptical NFW profiles fitted to the Al 15 weak-lensing data. Here and in Table 6, X and y positions are given with 
respect to the cluster catalogue positions in Table 1, with positive x corresponding to positive on the sky, i.e. negative in RA. 
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Figure 6. A1914. The maps at the top-left and top-right are from the SA before and after source-subtraction; in this and all subsequent SA figures, contours 
are every lo~, and follow levels —12, —10, —8, —7, —6, —5, —4, —3, 3, 4, 5, 10, 15, 20, 30. Sources marked with a '+' were modelled, while those marked 
with a 'x' were directly subtracted. The '□' marks the position of the centre of the fitted isothermal /3-model. The parameters of the labelled sources can be 
found in Table 8. Subfigure (c) shows A1914 weak-lensing 2D posteriors; in this and all subsequent weak-lensing 2D posterior plots, x and y are in units of 
arcseconds, c is dimensionless and M is total mass (integrated to f2oo)- in h~ Mg. Subfigure (d) shows A1914 SZ 2D posteriors. In this and all subsequent SZ 
2D posterior plots, x and y are in units of arcseconds, r c in kpc, /g as in h~ [ , and total mass (integrated to rioo as described in Section 4) in ft _I M s . The 
bottom-right image shows a composite X-ray, mass and SZ image. The greyscale shows Chandra X-ray data (Govoni et al. 2004) smoothed by a Gaussian of 
radius 15 arcsec, with 30% contours. Thicker black contours show the LENSENT mass reconstruction in 20% levels of density. Thick grey contours show the 
SA SZ decrement in 3o" intervals, where the map noise o" = 1 15 jliy. The white cross indicates the position, orientation and axis ratio of the elliptical /3-model 
fitted to the SZ data. The black cross indicates the maximum-likelihood position, orientation and axis ratio of the elliptical NFW profile fitted to the lensing 
data. 
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5.2 A1914 

This cluster is a fairly complex merger with two distinct X-ray 
clumps which appear to be moving in opposite directions east- 
west. At the same time, the overall mass distribution is irregu- 
lar and elongated NE-SW. Govoni et al. (2004) find filaments of 
hot gas connecting the mass peaks. Even before source-subtraction 
is performed, it is clear that A1914 has an elliptical shape in the 
SZ (Fig. 6a). It is surrounded by a large number of point sources 
which are clearly visible in the SA data. After source-subtraction 
(Fig. 6b), the majority of source contamination is removed cleanly. 
There are some 2o"sa residuals but these are more likely to be noise 
or faint extended sources than subtraction artifacts. The only signif- 
icant residuals are near two sources of flux< 4(Tla- The image here 
is consistent with that made during the early stages of science ob- 
servations with AMI (AMI Collaboration et al. 2006). 

Fitting an elliptical j3 -model to the data gives the parameters 
shown in Table 6; a circular model was also fitted but gave results 
with lower evidence. The centroid position is slightly offset south 
and east from the catalogue position, which agrees well with the 
overall structure of the X-ray data. It is also clear that SZ data 
traces out a much larger region of the extension of the gas, per- 
haps produced from the merger of the two sub-clumps visible in 
X-ray observations (Govoni et al. 2004). However, the individual 
clumps are not resolved. 

Lensing analyses of A1914 are difficult due to its highly 
distorted mass distribution. The contours of Fig. 6e show the 
LensEnt mass reconstruction from our optical data. The distribu- 
tion is clearly distorted and elongated, and looks similar to the dis- 
tribution reconstructed by Okabe & Umetsu (2008), but with lower 
resolution, so that the two different mass peaks remain unresolved. 

An elliptical NFW profile was fitted to the data, resulting in 
a mass estimate of 4.9 ± 1.4 x 10 14 /j _1 Mq, which is just in agree- 
ment with the lower end of Okabe & Umetsu's stated estimate of 
MvMai = 6.14 ±3.19 x 10 14 /i- 1 M Q , but we find that the mass is 
rather less concentrated than their estimate: c = 2.0 ± 1.3, as op- 
posed to c = 4.13 ± 2.79. These results are still consistent within 
the errors. The model that agrees best with the data has an elliptical 
geometry of 7] = 0.7±0.1 inclined at (j) = 150±30°. This fits with 
the mass reconstruction; essentially we are fitting a single profile to 
the distribution that Okabe & Umetsu separate into two peaks. 

An attempt was made to fit two components to this data; how- 
ever the evidence for two components was lower than that for a 
single ellipsoid. This could be due to a lower-resolution catalogue 
than that available to Okabe & Umetsu, who do see two mass peaks. 
However it also appears that their detected peaks overlie a mass el- 
lipsoid, rather than being two distinct clumps, so an ellipsoid might 
simply be a better description of the (clearly irregular) mass distri- 
bution. 

A reconstruction of the cluster mass and gas is shown in 
Fig. 6e along with the SZ data and the locations of the X-ray 
clumps. It is interesting to note from the SZ data that on the largest 
scales, aside from the high ellipticity, the gas is fairly uniformly 
distributed. Unlike the complex temperature-dependent structures 
picked out by the X-ray data, the SZ map provides a smooth and 
large-scale picture of the gas distribution. Of note is the north-west 
extension of the gas that is invisible in the X-ray data. 

In conclusion, the distributions of both gas and mass are 
highly irregular and detailed hydrodynamic simulations are nec- 
essary to improve understanding of the merger in A1914. 



5.3 A2111 

A2111 is widely agreed to be a head-on merger between two 
smaller clusters. The X-ray data presented by Wang et al. (1997) 
show an elongated morphology with two components. The subclus- 
ter was determined to have entered from the north-west, its core 
heating as it entered the gas in the centre of the main cluster. The 
outskirts of the cluster gas are fairly relaxed and the most disturbed 
gas lies only in the centre. 

Despite the large number of contaminating radio sources, 
the SZ decrement of A2111 is clearly visible in the source- 
unsubtracted SA map (Fig. 7a). A circularly- symmetric /3-model 
was fitted to the gas distribution using the SZ data, resulting in a 
centroid with parameters given in Table 6. 

A NFW profile was fitted to the lensing data: a circularly- 
symmetric model was slightly preferred over an elliptical model 
by the evidence values: the parameters from the resulting model 
are given in Table 6. 

Fig. 7e shows a composite image of the X-ray, lensing and SZ 
data. The gas and mass appear relaxed and circularly symmetric; 
the total mass is found to be 6.9± 1.1 x 1O 14 /i _1 M0 from the lens- 
ing data; this is consistent within the errors with the value measured 
via the SZ, 6.3±2.1 x 1O 14 /!~'M . Unlike X-ray observations and 
the long-baseline SZ observations by LaRoque et al. (2006), we are 
mapping the gas on the edges of the cluster, as well as the denser 
core. We thus find higher values for j3 and r c than are typical in the 
literature. 

5.4 A2259 

The SZ decrement of A2259 is clearly visible; indeed it is more 
striking than any of the nearby point sources (Fig. 8a). After these 
are fitted and removed, the decrement (Fig. 8b) appears slightly less 
elliptical and more reminiscent of the X-ray image: a composite is 
shown in Fig. 8d. The results of the /3 -model fitted to the SZ data 
are shown in Table 6. 

A2259 lies in the Galactic plane, and this was immediately no- 
ticeable when analysing the lensing data. Despite attempting vari- 
ous spatial, magnitude and colour cuts to the background galaxies, 
the evidence for a lensing effect around A2259 was very low, so that 
when fitting a NFW profile, the results from the posterior highly re- 
semble the priors. It appears that contamination from Galactic fore- 
grounds prevented us from detecting the small lensing signal from 
this low-mass cluster. 

5.5 A611 

A61 1 is the most dynamically-relaxed cluster in the sample, having 
a very uniform X-ray map with little substructure (LaRoque et al. 
2006). The 15.7-GHz source environment is clean; the decrement 
is immediately obvious even in the source-unsubtracted SA map 
(Fig. 9a). Fitting a j3 -model to the gas distribution results in the 
parameters in Table 6. 

Fitting a circularly- symmetric NFW profile to the lensing data 
results in the parameters in Table 6. An elliptical profile was also 
fitted, but the evidence was not higher than that for the circular pro- 
file, so the extra model parameters were not justified. The SZ and 
lensing masses agree well with each other at M ~ 5 x 10 14 /z~'Mq. 
Fig. 9e is a composite image of the X-ray data, LensEnt mass 
contours and SZ decrement. Okabe et al. (2010a) find a more dis- 
torted mass-distribution for this cluster but do see their mass distri- 
bution peak close to the position we find. 
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Figure 7. A2111. For the maps and posterior distributions, annotations and axis ranges and labels are as in Fig. 6. The parameters of the labelled sources in 
the SA maps can be found in Table 9. For the compositie image, the greyscale shows Chandra X-ray data smoothed by a Gaussian of radius 15 arcsec, with 
30% contours. Thicker black contours show the LENSENT mass reconstruction in 20% levels of density. Thick grey contours show the SZ decrement in 2o" 
intervals, where the map noise cr = 90,uTy. The white cross indicates the fitted centre of the circular /3-model with arm lengths of 1 <7 error on the position. 
The black cross indicates the fitted centre of the circular NFW profile with arm lengths of 1 cr error on the position. 
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Cluster 


A1914 


A2111 


A2259 


A611 


A851 


Parameter 


Lensing SZ 


Lensing SZ 


SZ 


Lensing SZ 


Lensing 
CI C2 


SZ 

CI C2 


x 1 arcsec 
y 1 arcsec 

c 

Ml lO^/r'Mg 
rj /i~'kpc 
P 


28 ± 36 43 ± 8 
-37±27 -17±7 

150±S 38^ 6 
0.7 ±0.1 0.7 ±0.1 

2.8 ±1.4 

4.9 ±1.4 5.5±1.5 

240 ±110 
1.4 ±0.4 


80±9 38±16 
68±11 50±16 

1.9 ±0.5 

6.9±1.1 6.3±2.1 
565 ± 195 
1.7 ±0.5 


38±21 
2±17 

400 ±160 
1.6 ±0.4 


47±9 30±15 
13 ±5 15 ±16 

4.0 ±1.3 

4.7 ±1.2 6.0 ±1.9 
490 ±230 
1.6 ±0.5 


33 ±5 107±5 
-18±2 -52±2 

5.1 ± 1.3 5.1 ±0.9 
3.2±1.0 3.4±1.0 


-88±5 115±30 
-27 ±4 120 ±25 

3.9±0.5 4.1 ±1.0 
155 ± 55 440+gg 
2.0±0.3 1.8±0.3 



Table 6. Mean posterior values for the j3-profile fits to the SZ data and NFW-profile fits to the lensing data for A1914, A2111, A2259, A611 and A851. For 
the last, values for the two different components fitted to each dataset are shown. 



Romano et al. (2010) also perform a weak-lensing analysis 
of A611 using data from the Large Binocular Telescope and we 
note that they also find a relaxed shape to the mass distribution. 
Fitting a NFW profile, they estimate M 20 o = 5.6+^ x 1O 14 M (= 
3.9+jJ x 10 14 fc _1 M©), with a concentration parameter c = 3.9^11 
and T2oo = 1545^3Qgkpc. Our measured mass and concentration 
parameter are higher but not significantly so given the errors (Ta- 
ble 6). We derive, straightforwardly from M200, ^200 = 1560 ± 
160 kpc, which is in extremely good agreement with Romano et al. 
despite the strong degeneracy between c and 7-200 in their model. 
Our measurement is also in agreement with Okabe et al. (2010a), 
who find M 2 oo = 5.47l"}-j] x 10 14 /! _1 M Q . These recent results and 
our own are lower than the estimate of Allen et al. (2003), who 
use Chandra data to calculate M200 = 6.6+5 3 j 3 x 10 14 /i~ 2 M Q (= 
9.4+^ 3 9 2 x lO'^-'Mg). 

5.6 A851 

A85 1 is the highest-redshift cluster in the sample, and is known to 
have a very irregular gas distribution, indicating that it is dynami- 
cally young: De Filippis et al. (2003) find that A851 is composed 
of two subclusters in the process of merging, the angular separa- 
tion of which is 1.5 arcmin. The SA map before source subtraction 
(Fig. 10a) shows a busy source environment, with several sources 
directly over the cluster. Despite this, the edges of the SZ decre- 
ment are still visible to the east and north-west of the largest con- 
glomeration of sources. Fitting 11 of the 15 sources (Table 12) and 
subtracting them from the SA data results in the source-subtracted 
map shown in Fig. 10b. The remaining unsubtracted sources are 
of flux densities comparable to the SA map noise. The SZ decre- 
ment is clearly visible and there are no significant positive residuals 
remaining on the map. The gas distribution is fairly distorted and 
ellipsoidal, with a possible extended spur to the north-west. 

A jS -model was fitted to these data: an elliptical model was 
preferred over a spherical model, but the model with the highest 
evidence was two spherical /3 -models at two separate positions: the 
resulting parameters from these fits are listed in Table 6. and the 2D 
posteriors are shown in Fig. 10c. 

The z = 2 cluster candidate detected by Dressier et al. (1993) 
might produce some SZ signal. Unfortunately its distance from the 
centre of the A85 1 decrement is less than 1 arcmin, so this cannot 
be resolved by the SA. 

Using LensEnt to reconstruct the mass distribution of A851, 
at a higher resolution than shown in Fig. 10, the two most signif- 
icant features are two ellipsoidal subclumps separated by around 
1 arcmin. These agree fairly well with the positions of the X-ray 
subclumps detected by De Filippis et al. They were also detected 



as mass overdensities in lensing analysis carried out by lye et al. 
(2000). 

Since two distinct components are clearly identified,we fit- 
ted two circularly-symmetric NFW profiles to the data. Elliptical 
models were also attempted but the evidence did not justify the 
extra parameters. Table 6 shows the resulting parameters, which 
indicate that the two subclusters are of roughly equal mass and 
concentration, and Fig. lOd shows the 2D lensing posteriors. The 
summed mass of these clumps, M = 6.6 ± 1.4 x lO 14 /i _1 M0, is 
consistent with the mass limits determined by Seitz et al. (1996): 
2.6 x 10 14 /j- 1 M Q < M < 7.1 x 1O 14 A _1 M0. Iye et al. (2000) do 
not attempt to make an estimate of the total mass of the cluster. 
Given that the redshift of this cluster is so high, we could be in- 
cluding more foreground (unlensed) galaxies in the field galaxy 
selection, reducing the signal and thus the detected mass, which 
could explain the 2a difference between our SZ and lensing mass 
measurements for this cluster. 

Fig. lOe shows a composite image of the SZ decrement, mass 
reconstruction from lensing and X-ray data. The north-east mass 
subclump appears coincident with the X-ray subclump, while the 
south-west mass concentration is a little south of the other X-ray 
subclump. The main body of the SZ decrement covers the X-ray- 
bright area and the two mass clumps, but there is also an exten- 
sion to the south-east, which is clearly visible even before source- 
subtraction. This might imply that gas is being 'squirted' out of the 
sides of the cluster, in directions perpendicular to that of the mo- 
tion of the subclumps. De Filippis et al. find that the compressed 
gas which also follows this NW-SE extension has a higher temper- 
ature of 6 keV; if the high temperature extends to the outskirts of 
the gas, this will also boost the observed SZ signal. 



6 DISCUSSION 

6.1 Mass comparison and modelling 

Fig. 1 1 shows a comparison of the different measured values of M 
for the clusters for which we were able to perform both the SZ and 
lensing analyses. Measured masses for A61 1 and A1914 are in very 
good agreement. Interestingly, we find better agreement between 
mass measurements for A1914 than Okabe et al. (2010b), who find 
that including this cluster significantly distorts their otherwise well- 
correlated M-T relationship. As we use SZ measurements at f^oo 
rather than X-ray at 7-500, this ma Y t> e indicative of SZ's advantage 
in measuring the mass of merging clusters, which are likely to be 
more common at the higher redshifts now being probed by new 
galaxy cluster surveys. 

We measure the mass of A851 as larger in the SZ compared 
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Figure 8. A2259. For the maps and posterior distributions, annotations and axis ranges and labels are as in Fig. 6. The parameters of the labelled sources in 
the SA maps can be found in Table 10. In the composite image, the grey scale shows Chandra X-ray data smoothed by a Gaussian of radius 15 arcsec, with 
30% contours. Thick grey contours show the SZ decrement in lo" intervals, where the map noise o" = 100 £iJy. The black cross indicates the fitted centre of 
the circular /3-model with arm lengths of lo" error on the position. 



to the lensing, although the difference is not a significant deviation 
given the errors. However it is the highest redshift cluster of the 
sample, it is possible that we are underestimating its mass via lens- 
ing, due to poor field galaxy selection. This may put limits on our 
ability to measure weak lensing at higher redshifts, at least with the 
two optical bands used in this analysis. 

For A21 1 1, the temperature may be somewhat higher than our 



M-T relation predicts, as this cluster is an ongoing merger and has 
a hot 8 keV central component; difficulties in reconciling disturbed 
clusters with theoretical relations have been seen in other studies 
(e.g. Cypriano et al. 2004). Given the small sample discussed here, 
it is difficult to draw general conclusions, but the technique would 
scale well to larger samples, from which it may be possible to cor- 
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Figure 9. A6 1 1 . For the maps and posterior distributions, annotations and axis ranges and labels are as in Fig. 6. The parameters of the labelled sources in the 
S A maps can be found in Table 1 1 . In the composite image, the greyscale shows Chandra X-ray data (LaRoque et al. 2006) smoothed by a Gaussian of radius 
15 arcsec, with 30% contours. Thicker black contours show the LENSENT mass reconstruction in 20% levels of density. Thick grey contours show the AMI 
SZ decrement in 2o~ intervals, where the map noise a = 80^Jy. The black cross indicates the cluster catalogue position. The white cross indicates the fitted 
centre of the circular /5-model with arm lengths of lo" error on the position. The black cross indicates the fitted centre of the circular NFW profile with arm 
lengths of 1 (j error on the position. 
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Figure 10. A851. For the maps and posterior distributions, annotations and axis ranges and labels are as in Fig. 6. Posterior distributions for two different 
components fit to each dataset are shown, indicated by the subscripts on the parameter values. The parameters of the labelled sources in the SA maps can 
be found in Table 12. In the composite image, the greyscale shows X-ray Multi-Mirror Mission -, Newton X-ray data, smoothed by a Gaussian of 5 arcsec, 
with 20% contours. Thicker black contours are 20% levels of a LENSENT reconstruction using a resolution of 100 arcsec. Thick grey contours show the SZ 
decrement in Iff intervals, where the map noise a = 75 Jijy. The black crosses mark the positions of the two isothermal beta profiles fitted to the SZ data, and 
the white crosses show the positions of the two NFW profiles fitted to the lensing data; the arm lengths show 1 o~ error bars in each case. 
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relate the state of merger with discrepancy from this M-T relation. 



6.2 Gas distribution on large scales 

The investigation of Navarro et al. (1995) into the suitability of the 
jS -model for X-ray clusters showed that the measured value of /3 
increases with increasing radius of fit. For relaxed clusters, typical 
values of j3 from X-ray and high-resolution SZ, both of which re- 
solve out large angular scales, are in the range 0.5 < j3 < 0.8. In the 
context of these core gas measurements, larger values of j8 usually 
indicate merging processes. 

LaRoque et al. (2006) provide the most recent combined SZ- 
X-ray analysis that covers the most relaxed clusters in our sample, 
A2111, A2259 and A611. We detect the gas masses in the similar 
proportions to LaRoque et al.: A61 1 and A21 1 1 are of similar mass, 
and more massive than A2259. This is reassuring, especially given 
that the hot (~ 8 keV) gas in the centre of A21 1 1 was ignored for 
our isothermal analysis. It is difficult to compare our values for r c 
and /3 since we are looking at a much more extended areas of the 
gas and one would expect the profiles to be different for different 
fitting areas. 

Using a /3 -model presents a problem if extrapolated to high 
radii: the density does not steepen quickly enough and large, non- 
physical gas masses are predicted. Therefore both high- and low- 
angular scale measurements specify a cut-off radius at which to 
measure the gas mass, often rgoo, while in this analysis we use j^oo- 
It can be difficult to extrapolate models produced from our data to 
smaller radii such as f25oo(~ 0.3r2oo), since with a high value of 
/3, our profiles may not actually reach that density before the mass 
goes to zero. 



6.3 Gas fraction 

For the most relaxed cluster in our sample, A611, we further in- 
vestigated the gas fraction using a joint SZ and lensing analysis. 
We removed the WMAP prior on gas fraction, giving it instead a 
uniform prior between 0.01 and Ih . The posterior plots for the 
model with and without the lensing data are shown in Fig. 12. To- 
gether, the datasets constrain the gas fraction, which cannot be done 



E 

B 5 



with either individually. It is noticeable that the SZ data has no bias 
in its measurement of the M — / gas degeneracy; the combined data 
posterior lies along the posterior from the SZ data alone. This im- 
plies our parameterization is good, and that with a sensible prior on 
/g as , SZ data allow us to make good estimates of cluster masses. 

The combined data favour a higher value of gas fraction for 
this cluster than the WMAP prior normally used: 0.23 ±0. \h 
compared to 0. 123 ± 0.02h~ l . However the errors are large so it is 
difficult to tell whether this is significant. It would be illuminating 
to extend this analysis to a large dataset of relaxed clusters, as gen- 
eral conclusions about cluster morphology cannot be made from 
this single result, and the other lensing-detected clusters in this pa- 
per are complex mergers. 



6.4 Data issues and cluster selection 

This pilot study has been helpful in identifying some areas of dif- 
ficulty in joint SZ-weak-gravitational-lensing analysis. The high 
density of foreground optical sources in the field of A2259 made 
measurement of its low mass difficult via gravitational lensing, 
while it was well-measured via the SZ. On the other hand, the SZ 
analysis of A115 suffered from the effects of contaminating radio 
sources, and this effect is also seen in other AMI cluster observa- 
tions, e.g. by AMI Consortium: Shimwell et al. (201 1). 

For those clusters observed outside the Galactic plane, with 
unresolved (easily-modelled) radio sources, the analysis showed re- 
markable agreement between measured M200 regardless of dynam- 
ical state, X-ray luminosity, temperature, redshift and number of 
detected sub-structures. Currently the AMI Consortium has over 50 
further detections of SZ clusters in similar environments covering a 
redshift range of 0.0894 ^ z ^ 0.686 and an X-ray luminosity range 
of 3.0 x 10 37 < L x ^ 28 x 10 37 W, from X-ray-selected samples. 
Of these, 11 have associated archival CFHT data, on which may 
be performed a similar analysis. Also, the clusters being discov- 
ered by Planck (Planck Collaboration 201 la) are highly amenable 
to observations by AMI (AMI Consortium: Hurley- Walker et al. 
2011), so there is potential to expand the sample size by another 
order of magnitude over coming years. 



7 CONCLUSIONS 

(i) Using AMI SZ data that measure out to r2oo, an d thus probe 
scales comparable to those of weak-lensing observations, we have 
used a fast, Bayesian analysis to produce posteriors for useful pa- 
rameters for six clusters, on the assumptions of an isothermal j8- 
model to SZ data and an NFW model for the lensing data; 





i 8 10 12 14 16 18 
M/IT 1 M ra ,,„« 



Figure 11. Comparison of the different values of M obtained by the lensing 
and SZ analyses. The dashed line shows a 1:1 ratio. 



Figure 12. / gas -M degeneracy for SZ data alone (left) and SZ & GL data 
combined (right). 
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(ii) Of the four clusters for which we have both weak lensing 
and SZ data, we find that the mass estimates are in very good agree- 
ment, and that we may need to improve our field-galaxy selection 
for clusters at z > 0.4 for a larger sample; 

(iii) We perform the first multiple-component weak-lensing 
analysis of A115 and discover significant substructure; 

(iv) We confirm the unusual separation between the gas and 
mass components in A1914; 

(v) For A611, the most relaxed cluster in the sample, we have 
carried out a joint weak-lensing and SZ analysis with the gas frac- 
tion as a free parameter, and find that / gas = 0.23 ± 0. lO/i -1 . 
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ID 


RA 


Dec 


5 L A/mJy 


o:la 


ID 


RA 


Dec 


5 LA /mJy 


ccla 


A 


00 56 55.18 


+26 31 24.48 


52.0±2.8 


-0.5 ±0.3 





00 56 01.38 


+26 28 43.78 


1.7 ±0.3 


1.1 ±1.4 


B 


00 56 02.85 


+26 27 20.98 


27. 6± 1.4 


0.8±0.1 


P 


00 56 07.16 


+26 29 46.05 


1.4 ±0.3 


0.5± 1.5 


C 


00 57 21.47 


+26 17 27.90 


7.1 ±0.7 


0.9± 1.1 


Q 


00 56 06.97 


+26 30 51.52 


1.3 ±0.3 


0.4± 1.8 


D 


00 55 50.46 


+26 24 39.02 


5.4 ±0.3 


2.5 ±0.4 


R 


00 56 00.94 


+26 25 03.79 


1.1 ±0.2 


1.3 ± 1.3 


E 


00 56 20.93 


+26 30 55.24 


4.1 ±0.3 


1.1 ± 1.0 


S 


00 55 59.02 


+26 17 49.46 


1.0±0.1 


1.3 ± 1.3 


F 


dn n£ 
UU J / Oo. j j 


i 15 fY2 
+ZO ZJ Z/.UJ 


4.0 ± 0.5 


1.5 ± 1.4 


T 


t\r\ 1 f\ CA 
UU JO 1U.04 


n nl /in 
+ZO Iz UJ.4U 


1 .0 ± 0. 1 


0.5 ± 1 .4 


G 


00 56 08.49 


+26 25 13.46 


3.6 ±0.2 


-0.5 ±0.6 


U 


00 56 28.46 


+26 24 03.69 


0.9 ±0.1 


1.3 ± 1.5 


H 


00 56 22.47 


+26 23 03.43 


2.8 ±0.2 


0.5±0.6 


V 


00 55 58.37 


+26 25 33.95 


0.9 ±0.2 


0.8±1.7 


I 


00 56 55.38 


+26 29 12.30 


2.7 ±0.6 


0.1 ±1.6 


w 


00 55 48.04 


+26 25 56.77 


0.7 ±0.1 


2.3± 1.5 


J 


00 56 18.52 


+26 05 00.41 


2.3 ±0.3 


-0.9 ±1.3 


X 


00 56 04.35 


+26 25 00.45 


0.7 ±0.2 


-0.9±1.5 


K 


00 54 55.80 


+26 19 13.30 


2.2 ±0.3 


-0.3 ±1.4 


Y 


00 56 05.32 


+26 23 44.74 


0.6 ±0.1 


0.1 ± 1.6 


L 


00 56 50.45 


+26 29 38.97 


2.1 ±0.5 


0.4±1.6 


Z 


00 56 10.69 


+26 24 15.30 


0.6 ±0.1 


2.2±1.5 


M 


00 56 43.35 


+26 29 22.20 


2.0 ±0.4 


0.2±1.5 


AA 


00 55 49.68 


+26 22 20.83 


0.5 ±0.1 


0.2±1.6 


N 


00 56 40.69 


+26 30 10.21 


1.8 ±0.4 


2.2±1.5 













Table 7. Sources found in the LA image of A115, in order of LA flux density. In this and all subsequent tables, the flux densities given are measured from 
the combined-channel 15.7-GHz LA maps. The spectral index is calculated by measuring the flux densities from the separate channels and fitting a power-law 
spectrum, lo" errors are shown for both measurements. 



ID 


RA 


Dec 


5 L A/mJy 


o:la 


SsA/mJy 


asA 


A 


14 27 25.06 


+37 46 33.19 


10.7 ±0.5 


0.1±0.6 


8.3 ±0.7 


0.2±0.5 


B 


14 25 08.41 


+37 52 42.58 


4.3±0.2 


1.7±0.7 


4.4 ±0.3 


1.3±0.5 


C 


14 25 52.42 


+38 03 04.66 


2.9±0.3 


-0.5 ±1.2 


3.8 ±0.3 


-0.1 ±1.0 


D 


14 25 40.83 


+37 45 47.93 


2.7±0.1 


1.2 ±0.5 


3.1 ±0.2 


1.0±0.4 


E 


14 25 05.12 


+37 55 15.84 


2.7±0.3 


0.8±1.3 


2.5 ±0.3 


1.0± 1.1 


F 


14 27 10.92 


+37 55 13.89 


2.7±0.3 


0.9±1.3 


2.5 ±0.4 


1.6± 1.2 


G 


14 25 57.37 


+38 01 10.48 


1.4±0.3 


0.3±1.6 


1.2±0.2 


0.2 ±1.4 


H 


14 25 47.75 


+37 47 48.32 


1.2±0.1 


-0.4 ±1.2 


1.2±0.1 


0.0±1.0 


I 


14 25 56.73 


+37 55 10.47 


0.8±0.1 


0.7±1.7 


0.5 ±0.1 


0.1 ±1.4 


J 


14 25 50.36 


+37 45 07.81 


0.7 ±0.1 


0.8±1.5 


0.7 ±0.1 


1.1 ±1.5 


K 


14 25 05.78 


+37 47 02.26 


1.2±0.2 


-0.5 ±1.7 






L 


14 25 43.15 


+37 40 06.15 


0.8±0.1 


-0.8±1.7 






M 


14 25 58.12 


+37 43 58.58 


0.6±0.1 


1.0=b 1.5 






N 


14 25 38.95 


+37 57 32.50 


0.6±0.1 


1.6±1.7 







Table 8. Sources found in the LA image of A1914, and fitted to the SA data, in order of LA flux density. In this and all subsequent tables, the fitted SA flux 
densities and spectral indices with lo~ errors are shown in the final two columns; sources without entries in these columns were directly subtracted from the 
data to reduce the parameter space to a more manageable size. 
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ID 


RA 


Dec 


SLA/mJy 


aLA 


5 SA /mJy 


a S A 


A 


15 40 55.25 


+34 30 17.50 


17.4±0.7 


-0.7 ±0.4 


14.6±0.5 


-0.1 ±0.2 


B 


15 38 46.64 


+34 18 58.04 


6.3±0.3 


1.6 ±0.6 


5.5 ±0.3 


1.5 ±0.5 


C 


15 40 31.18 


+34 30 09.72 


2.9±0.2 


0.2± 1.1 


2.9 ±0.2 


1.0 ±0.7 


D 


1 c m 1 1 QH 

1 J jy 11.5/ 


i 1A on 11 fYl 

+J4 ZV dj.V/ 


2.5 ±0.1 


0.7 ± 0.7 


i /; _i_ n 1 
2.0 ± 0. 1 


1 c _l n a 
I.J ± 0.4 


E 


15 39 55.06 


+34 20 11.40 


2.3±0.1 


1.2 ±0.6 


1.9±0.1 


1.2 ±0.5 


F 


15 40 41.89 


+34 18 35.72 


2.1 ±0.4 


-0.3 ±1.3 


2.8 ±0.2 


0.7 ±1.0 


G 


15 40 49.55 


+34 30 35.66 


1.9±0.6 


0.8 ±1.6 


1.9 ±0.3 


0.7 ±1.5 


H 


15 40 31.90 


+34 28 16.52 


1.5±0.3 


0.2 ±1.4 


2.3 ±0.2 


0.3 ±1.0 


I 


15 38 49.51 


+34 26 56.26 


1.3±0.2 


1.7± 1.5 


1.8±0.2 


2.3±1.1 


J 


15 39 08.11 


+34 21 09.02 


0.8±0.1 


0.4 ±1.6 


0.7 ±0.1 


0.2 ±1.4 


K 


15 39 56.78 


+34 29 33.38 


0.8±0.1 


0.8 ±1.4 


0.9 ±0.1 


-0.5 ±1.0 


L 


15 39 30.11 


+34 29 05.46 


0.6±0.1 


1.2±1.6 


0.5 ±0.1 


1.0±1.4 



Table 9. Sources found in the LA image of A2 1 1 1 , and fitted to the S A data, in order of LA flux density. 



ID 


RA 


Dec 


S LA /mJy 


«LA 


S SA /mJy 


«SA 


A 


17 19 13.69 


+27 44 49.47 


2.2±0.3 


2.3±1.3 


1.9±0.3 


2.4±1.1 


B 


17 19 30.18 


+27 46 16.12 


1.1 ±0.2 


0.5 ±1.6 


0.9±0.2 


0.5 ±1.4 


C 


17 20 51.25 


+27 45 08.05 


1.4±0.1 


-0.2±1.5 


1.6±0.2 


0.7 ±1.2 


D 


17 19 27.06 


+27 32 51.05 


1.0±0.2 


0.1 ±1.7 


0.6±0.2 


0.0 ±1.6 


E 


17 20 20.26 


+27 28 33.96 


1.0±0.2 


1.3±1.8 


1.2±0.2 


2.0 ±1.5 


F 


17 20 07.54 


+27 44 35.03 


0.8±0.1 


0.1 ±1.3 


0.6±0.1 


0.2 ±1.2 


G 


17 19 44.57 


+27 32 37.99 


0.6±0.1 


1.7±1.6 


0.4±0.1 


1.4± 1.5 


H 


17 19 41.33 


+27 42 10.07 


0.4±0.1 


0.4 ±1.7 


0.4 ±0.1 


0.6±1.7 


I 


17 20 16.68 


+27 37 50.27 


0.3±0.1 


1.0± 1.7 


0.5 ±0.1 


0.4 ±1.4 



Table 10. Sources found in the LA image of A2259, and fitted to the SA data, in order of LA flux density. 



ID 


RA 


Dec 


5 L A/mJy 




SsA/mJy 


«SA 


A 


08 00 07.83 


+36 04 08.56 


5.9±0.2 


0.5 ±0.4 


6.1±0.2 


0.4±0.3 


B 


08 02 00.53 


+36 08 54.64 


3.3±0.3 


-0.2 ±1.0 


2.6 ±0.2 


0.3±0.9 


C 


08 00 40.55 


+36 14 23.78 


3.1±0.3 


0.3± 1.1 


3.6±0.3 


1.2±0.7 


D 


07 59 51.35 


+36 11 05.49 


3.0 ±0.4 


0.2±1.2 


2.5 ±0.4 


0.8± 1.1 


E 


08 00 43.25 


+36 14 03.44 


3.0±0.3 


0.3± 1.1 


2.8 ±0.3 


0.4±0.8 


F 


07 59 55.92 


+35 58 35.84 


2.8±0.2 


1.0± 1.3 


1.5 ±0.2 


0.6± 1.1 


G 


08 02 12.27 


+36 03 48.26 


2.2±0.2 


1.0± 1.3 


1.4 ±0.3 


1.2± 1.2 


H 


07 59 48.28 


+36 06 41.42 


2.0±0.2 


0.2±1.5 


2.2 ±0.3 


-0.4±1.3 


I 


08 00 03.25 


+36 00 51.74 


1.0±0.2 


-0.5 ±1.6 


0.2 ±0.1 


-0.3±1.6 


J 


08 00 59.26 


+35 55 51.49 


0.9±0.1 


0.3 ±1.3 


0.9 ±0.1 


0.6± 1.1 


K 


08 01 17.11 


+36 04 31.27 


0.6±0.1 


1.7± 1.4 


0.5 ±0.1 


0.7± 1.1 


L 


08 00 30.27 


+36 00 41.98 


0.6±0.1 


1.1 ±1.6 


0.7 ±0.1 


0.4 ±1.4 


M 


08 01 24.68 


+36 05 37.48 


0.6±0.1 


-0.0 ±1.5 


0.5 ±0.1 


-0.1 ± 1.3 


N 


08 00 52.72 


+36 06 13.52 


0.4±0.1 


0.7 ±1.8 


0.3 ±0.1 


0.7 ±1.6 


O 


08 00 40.14 


+35 59 51.55 


0.4±0.1 


1.3 ± 1.6 


0.3 ±0.1 


1.7± 1.5 


P 


08 00 11.88 


+35 50 15.31 


1.6±0.3 


-0.9±1.6 






Q 


08 00 54.95 


+36 17 43.76 


1.4 ±0.4 


-0.8±1.7 






R 


08 00 38.30 


+36 10 56.43 


0.6±0.1 


-0.6±1.6 







Table 11. Sources found in the LA image of A611, and fitted to the SA data, in order of LA flux density. 
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Bayesian analysis ofweak-gravitational-lensing and Sunyaev-Zel 'dovich data for six galaxy clusters 



ID 


RA 


Dec 


5 LA /mJy 


aLA 


SsA/mJy 




A 


09 44 48.36 


+47 00 02.10 


6.0±0.3 


0.4 ±0.6 


1.7 ±0.8 


0.0± 1.7 


B 


09 42 57.44 


+46 58 49.93 


2.2±0.1 


0.9 ±0.5 


2.9 ±0.5 


0.5±0.6 


C 


09 42 34.86 


+47 18 23.55 


2.1 ±0.5 


-0.1 ± 1.7 


2.2 ±0.1 


0.7 ±0.3 


D 


09 43 32.44 


+46 47 19.16 


1.4 ±0.3 


— 0.7 ± 1.5 


1.2 ±0.2 


-0.5 ± 1.3 


E 


09 43 19.30 


+46 47 32.98 


1.2±0.3 


0.1 ± 1.6 


1.6 ±0.2 


-0.1 ± 1.2 


F 


09 43 12.10 


+47 04 22.51 


1.2 ± 0.1 


1.9±0.8 


1.2±0.1 


0.9±0.8 


/ ■ 

Lj 


no /io O/i 1 n 
uy 4Z Z4. lu 


_l/1"7 m A Q 1 Q 
+*+ / UZ ^fo. Ly 


1 . 1 ± u.u 


n o 4- n o 
u.y ± u.y 




1 14.(1 ^ 
1 .Z It U.J 


H 


09 42 23.62 


+46 50 07.13 


1.0 ±0.2 


0.3 ±1.4 


1.4 ±0.2 


-0.2±1.2 


I 


09 42 06.14 


+46 57 05.32 


0.9±0.1 


0.9 ±1.7 


0.6 ±0.1 


1.1 ± 1.5 


J 


09 42 49.52 


+46 57 07.06 


0.6±0.1 


2.1 ±1.4 


0.3 ±0.1 


1.6± 1.2 


K 


09 43 00.70 


+47 01 27.71 


0.6 ±0.1 


1.1 ±1.4 


0.6 ±0.1 


1.5±1.0 


L 


09 44 05.15 


+46 44 02.43 


4.5 ±0.5 


0.3 ±1.4 






M 


09 42 19.30 


+46 56 55.84 


0.5±0.1 


-0.7 ±1.7 






N 


09 43 44.18 


+47 02 54.42 


0.4±0.1 


-0.3±1.7 






O 


09 42 31.73 


+46 55 07.53 


0.4±0.1 


-0.3 ±1.6 







Table 12. Sources found in the LA image of A851, and fitted to the SA data, in order of LA flux density. 
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